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Abstract 

We show that solutions to Smoluchowski's equation with a con- 
stant coagulation kernel and an initial datum with some regularity 
and exponentially decaying tail converge exponentially fast to a self- 
similar profile. This convergence holds in a weighted Sobolev norm 
which implies the L 2 convergence of derivatives up to a certain or- 
der k depending on the regularity of the initial condition. We prove 
these results through the study of the linearized coagulation equation 
in self-similar variables, for which we show a spectral gap in a scale of 
weighted Sobolev spaces. We also take advantage of the fact that the 
Laplace or Fourier transforms of this equation can be explicitly solved 
in this case. 
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1 Introduction 

1.1 Smoluchowski's equation with a constant 
kernel 

Smoluchowski's equation is a well-known model for the time evolution 
of irreversible aggregation processes [31 [2j [131 tH]- The object of our 
study is the continuous version of this equation, with the additional 
assumption that the rate at which any two nuclei coalesce (the coag- 
ulation rate or kernel) is independent of their sizes. If the number 
of nuclei of size y > at a given time t is given by f(t,y)dy, then 
Smoluchowski's equation reads: 

d t f(t,y) = C(f(t, •),/(*, -))(y), (1) 

where 

1 rv r°° 
C(f, f)(y) ■= = / f(x)f(y -x)dx- f(y) / f{x) dx. (2) 
z Jo Jo 

When other coagulation rates are considered, a coagulation coefficient 
a(x, y) appears in the equation; here we have set a(x, y) = 1 for all 
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x,y > 0. The integral J Q °° f(t,y)ydy, called the mass of f at time t, 
represents the total mass of all nuclei at time t, and is a conserved 
quantity for this equation. 

One of the more remarkable properties of Smoluchowski's equation 
is the self-similar behavior that its solutions exhibit at large times, 
which is expected by formal arguments but has been rigorously proved 
only in very particular cases. This behavior is our main concern in 
this paper, and among our main results we highlight the following one: 

Theorem 1.1. Take a locally absolutely continuous function /o : 
(0, oo) — ► [0, oo) such that, for some v > 0, 

POO 

/ Cf^))Ve^y<oc. 
J o 

Let f be the solution of Smoluchowski's equation (H|) with initial con- 
dition /o (this problem is known to be well-posed; see \lb\ \14V- Then 
there is a number p > such that for any < 5 < 1, 

POO 

/ (t 2 d y (f(t,ty))-d y g p y y A eM dy < K 2 r 25 forallt>0, (3) 
J o 

for some K > which depends only on the initial condition /o. In 
particular, 

sup {y 2 \t 2 f(t, ty) - g p (y) \ } < Kt' 5 for all t > 0, 
y>0 

— — 

where p is the mass of fo and g p {y) := -e p y for y > 0. 

This theorem is a direct consequence of Theorem 11.6 1 below . proved 
in Section [H Results along these lines were already obtained in [H] 
and [Td] , where it was proved that with some regularity of the initial 
condition it holds that 

sup{y \t 2 f(t,ty) - g p {y)\) -> when t -» oo. (4) 

y>0 

A result similar to the latter, but giving uniform convergence only 
in compact sets, can be found in [9]; weak convergence was proved 
in [12] using entropy arguments. Hence, the main improvement of 
Theorem 11.11 is that we give an explicit rate of convergence, which is 
to our knowledge a new result, and that the convergence in eq. (|3|) 
is stronger than uniform convergence. In fact, we give corresponding 
results for the rate of convergence of higher derivatives of / as long 
as the initial condition is regular enough. We note that this rate is 
optimal in general, as can be seen from the evolution of the moment 
of order of /, J f. These results, and the methods we use to arrive 
at them, are best described in terms of a rescaled form of eq. (pQ), 
which we introduce next. 
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1.2 The self-similar equation 



A scaling solution or self similar solution of Smoluchowski's equation 
is a solution / to eq. ([T]) which is a rescaling of some fixed function g 
for any time t: 

f(t,y)=p(t)g(q(t)y) (t,y >0), (5) 

where p, q > are functions of time. When such a solution exists, g 
is called a scaling profile or self-similar profile. 

Let us briefly describe the problem of the scaling behavior of 
Smoluchowski's equation for a general coagulation rate, and then we 
will focus on the problem for a constant rate. It is expected that for 
a general class of homogeneous coagulation kernels, there is a unique 
self-similar profile with given mass p > 0, and that for a very wide 
class of initial conditions solutions to Smoluchowski's equation approx- 
imate, for large times, the self-similar solution with the same mass, 
in a sense to be precised. For general homogeneous coagulation rates 
there are no rigorous proofs of this behavior, except for recent results 
that have shown the existence of self-similar profiles [6] and given 
some of their properties [5J 0]. Nevertheless, for the specific coagu- 
lation rates given by a(x,y) = x + y and a(x,y) = 1 (our case), it 
is known that the self-similar profile is unique for each given finite 
mass, and the convergence to it has been proved in [T31 HSl [12] - For 
a(x,y) = 1, the profile corresponding to a mass p is the function g p 
mentioned in Theorem ll.il and the convergence to it holds at least in 
the sense of eq. under quite general conditions on the initial data. 

In our case, if / is a self-similar solution of eq. (pQ) , then it is known 
that the scaling in eq. (|5|) is determined: it must happen that for some 
r G M, 

^ + ^ = (ir^iT7) for ^°^>°- 

Hence, up to time translations, we can consider that all self-similar 
solutions are of this form. The following change of variables is then 
useful: if / is a solution of eq. ([T]) for t G [0, oo), we set 

g(t, y) := e 2t f{e l - 1, e l y) for t > 0, y > 0, (6) 

whose inverse is 

f(t,y) = j^2 9^og(l + t), ^) fort>0,y>0. (7) 

The fact that / is a solution of eq. ([I]) is equivalent to g satisfying the 
following self-similar coagulation equation: 

d t g=2g + yd y g + C{g,g) (8) 
=:D(g) + C(g,g). (9) 
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Through this change of variables, self-similar solutions / of eq. ([T|) 
correspond, up to a time translation, to stationary solutions of eq. ([H]). 
Given this equivalence, we study eq. ([8]), which is more convenient 
when considering the scaling behavior of solutions. 

Let us be more precise as to what we mean by a solution to eq. (|SJ) . 
Below we denote our functional spaces by 

L 1 1=^(0, oo) L\ ■=!} ((0,oo); ydy). 

Definition 1.2. For any g,h £ L , we define C(g, h) by 

C(g,h)(y) :=~ [" g(y')h(y - y') dy dy' 
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g(y)h(y')dy' (10) 

h(y)g(y') dy'. 
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Notice that this is in agreement with the expression in eq. ©, and 
that this operator is symmetric: C(g,h) = C(h,g) for any g, h. 

Definition 1.3. A solution of eq. ([8]) on [0, oo) is a nonnegative 
function g G C([0, oo); L 1 PiL\) such that eq. (JS]) holds in the following 
sense: for all (ft € Cg°, # satisfies 



g(t,y)4>(y)dy = g(0,y)<j)(y)dy + / g{s,y)(j){y) dy ds 
i) Jo Jo Jo 

t POO 

/ yg(s,y)(p'(y)dyds 

JO 

+ / C(g(s),g(s))(y)(j)(y)dyds for t > 0. (11) 

JO JO 

The results in [3] or [H] prove that for any nonnegative g° G L 1 nL|[ 
there exists a unique solution to eq. © in the above sense with the 
initial condition g(0,y) = g°{y) for y > 0. 

Eq. ([8j) has a one-parameter family of stationary solutions {g p } p >o 
which is explicitly given by 

g P {y) :=-e~^ . (12) 

The solution g p has mass (first moment) p and number of particles 
(zeroth moment) equal to 2. If we write the result in eq. @ in terms 
of a solution g to eq. ([8]), then it means that 

lim sup{y|#(t,y) - g p (y)\} = 0. 

t^oo y>0 
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Some regularity conditions need to be imposed on the initial data for 
the latter convergence to hold; we refer the reader to [15], where the 
result is proved, for details. 

1.3 Main results 

Our main objective is to give an explicit exponential rate of conver- 
gence to equilibrium for solutions to the above equation in various 
norms. Explicitly, for an integrable function h : (0, oo) —* R, a posi- 
tive number p and an integer k > — 1 we consider the norms 

poo 

INlL : = / y 2(k+1) (D k h(y)) 2 e^dy, (13) 
J o 

where D k h represents the k-th derivative of h and D~ l h denotes the 
primitive of h given by 

/•oo 

D- l h{y) = -H{y):=- I h(x)dx. (14) 

Jy 

Note that the subindex p indicates the exponential weight, and not 
the power of the function as is usual when considering L p spaces. 
The latter norm with p = 2/p is directly suggested by the quadratic 
approximation to the following relative entropy functional near a sta- 
tionary solution, and is the initial motivation for it: 

F[g\g p ] := J™ (g{y) (log j&L - l) + g p (y)) <*V- (15) 

We prove both local and global convergence results; among the 
local ones we have the following, proved in Section [5j 

Theorem 1.4. Let g be a solution to eq. ([8]) with mass p and non- 
negative initial condition g° at t = 0. Take 0<fi<u<2/p and an 
integer k > — 1. Choose e > 0, and assume that 

1. The exponential moment of order v/2 is bounded by a constant 
E > for all t > 0; this is, 

roo 

I \g(t,y)\e% y dy < E forallt>0. 
Jo 

2. The relative entropy F\g®\g p \ is less than e. 

Then, for any < 5 < 1 there is an eo = eo(p,E,5,k, p,u) such that 
if e < eo then 

ll^-^II^^^H^-^ll^e- 5 * fort>0, (16) 
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for some constant K = K(p, E, 5, k, p, v). When k = —1, the same is 
true with \\g° — g P \\_ l ^ instead of \\g° — g P \\_ l v - 

Hs(*) " 9 P \L^ < K V ~ 9 P \\_ X ^ e' 5t for t > 0. 

Observe that these results are meaningful when \\g°\\ ku , or \\g°\\ k in 
this last case, are finite. 

Remark 1.5. Below we give intrinsic conditions on the initial data 
which are equivalent to point 1 in the above theorem (see Section 0J 
and concretely eq. ()92fl ). so it is really a condition on g°. In particular, 
it is satisfied if g (y) < ve~z y for all y > 0. 

Our main global convergence result can be stated as follows: 

Theorem 1.6. Let g be a solution to eq. (|8|) with mass p and non- 
negative initial condition g° such that (g )' € L . Assume that, for 
some v > and k > —1, 

/•oo 

/ 9°(y)e uy dy < oo and \\g°\\ < oo. 
Jo 

Then there exists < \i < 2/p such that for any < 5 < 1, 

\\g{t)-g p \\ k ^<Ke- St fort>0 (17) 
for some number K which depends in an explicit way on k, u, 5 and 

g°- 

For k = 1, this implies Theorem 11.11 through the self-similar change 
of variables in eq. (J6]). When k = —1 or 0, the condition ||5 ,0 || fc v < oo 
is redundant, as it is already implied by J g°e uy < oo and g' G L . 
We also remark that all the constants in the previous theorems are 
constructive and can be explicitly given by following their proofs. 

We point out that all the arguments in the paper can be in fact 
carried out with the norms 

/•oo 

INlL : = / y 2k (D k h(y)) 2 e^dy (18) 
Jo 

for k > and p > 0, and 

POO 

\ih\W\, = \\h\\\, := / H(yfe^dy, (19) 

J 

as before. Notice that here the power weight is lower in all norms 
with k > 0. In this case, Theorems 11.41 and 11.61 are still true with 
the sole difference that one obtains an exponential rate of convergence 
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with 5 < 1/2 instead of 5 < 1. In particular, one obtains a statement 
as that in Theorem 11.11 with a modified power weight and rate of 
convergence: in the conditions of Theorem ll.il and for < 5 < 1/2, 




(t 2 d y (f(t,ty))-d y g P ) y 2 e^dy 



< KH 



2 + -2<5 



for all t > 



and in particular 



sup{y\t 2 f(t,ty)-g p (y)\} < Kt' 6 for all * > 0. 



y>0 



We have not included the calculations for the norms (|18f) with different 
power weights in order to make the paper more readable, but there 
are no difficulties in rederiving our results in this case. 

We also note that the result on exponential convergence in the I? 
norm, with an exponential rate of convergence of 1/2, can be arrived 
at by explicit computations on the solution of the Fourier transform 
of eq. (H|). We do this in Lemma 16.11 

These results strongly suggest that the behavior near y = of 
solutions to eq. ([T]) is decisive when studying the speed of convergence 
of solutions to a self-similar one: it seems likely that one cannot obtain, 
for example, the convergence of derivatives at y = 0, so some weight 
near y = is probably unavoidable. 

The strategy of proof for our results stems from the relative entropy 
functional defined in eq. (|15p , which was first introduced in the context 
of eq. d8|) in [12], where it was proved that 



for any solution g of eq. (jSJ), where G and G p are the primitives 
of g and g p defined as in eq. (UH). That is: the above expressions 
are Liapunov functionals for this equation. If one wants to carry out 
the classical method of studying the convergence to equilibrium in 
a neighborhood of a stationary solution by looking at the spectral 
properties of the linearization of the equation near such a solution, 
the above functionals suggest that we consider the linearized operator 
in the norm which appears in their quadratic approximation. For eq. 
(l20l) . this norm is precisely ||-|| 2; for the functional in eq. OTj) . it is 



d F[g(t,-)\g p ] 
dt fg{t,y)dy 



< 



and 



(20) 




(21) 




2 . We follow this path and study the equation 



d t h = Lh := 2h + ydyh + 2C{h, g p ) 



(22) 
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where L is the linearization of the operator D(g) + C(g,g) from eq. 
([9|) near the stationary solution g p . We find that L, when restricted 
to functions h such that / °° yh(y) dy = 0, has a spectral gap in the 
norm \\-\\ k p for < fi < 2/p and integer k > — 1; see Sections 13.21 - 
13.41 for a precise statement (see also |5] for a spectral study of L in 
an L 1 framework). This shows exponential convergence to for the 
linear equation (|22p in these norms, and can also be used to prove 
exponential convergence to equilibrium for the nonlinear equation (JSJ) 
in a close-to-equilibrium regime. 

In order to treat solutions far from the equilibrium, the results of 
|15j were already available. However, as we are interested in the rate 
of convergence, we prove a similar result which gives an explicit rate 
of convergence in the L 2 norm, using estimates on the explicit solution 
for the Fourier transform of (JSj) which resemble those in the proof of 
uniform convergence in [15] . With further estimates on the derivatives 
of the explicit solution one could probably reproduce "by hand" our 
results on the rate of convergence for the derivatives. Nevertheless, the 
calculations soon become very involved, and in our opinion do not add 
much to the understanding of the equation, while a linear study proves 
to be less technical, and possibly better suited for generalization to 
other coagulation kernels. 

The paper is organized as follows: in Section [2] we derive our main 
functional inequality, obtained from a linearization of an inequality of 
Aizenman and Bak [lj, which is then used to prove that the linear 
operator L has a spectral gap in the norm ||-||_ 12 /p- I n Section [3] 
we define precisely the linear operator L in the spaces associated to 
the norms mentioned above, and prove that it has a spectral gap 
in them. Section 0] can be read independently from the rest of the 
paper; in it, we find an explicit expression for the evolution of the 
exponential moments of solutions to eq. ([H]). This is a remarkable 
piece of information, which can be obtained because of the widely 
known fact that the Laplace transform of eq. (H|) is explicitly solvable. 
In turn, this information on the exponential moments is needed in the 
arguments of later sections. In Section [5] we give local exponential 
convergence results for the nonlinear equation (JEj) with the help of the 
linear study from Section [3l and in particular we prove Theorem 11.41 
Finally, in Section [6] we state and prove our global convergence results, 
and in particular prove Theorem 11.61 

2 Functional inequalities 

The study of the linearization of eq. ([H]) in the norm 2 / p defined 
in eq. (fT3l) is directly suggested by the Liapunov functional in eq. ([2T]) . 
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Here we prove the necessary inequality to show that the linearized 
operator L (which we define more precisely later) has a spectral gap 
in the norm ||-||_ 1 2 / p \ m ore precisely 

(h,Lh)_ X 2 <-|HI_i2 (23) 

for any sufficiently regular function h with J °° yh(y) dy = 0, where 
(•>•)_! 2/p denotes the scalar product associated to the norm IHI-12/p' 
see Proposition 12.61 and Lemma 13.12 1 for a more precise statement. 

In order to prove this we use an inequality of Aizenman and Bak, 
which was first introduced in their paper [I] to study the convergence 
to equilibrium for eq. (pQ) with an additional fragmentation term. In 
|12j . Aizenman-Bak's inequality is also needed in the proof that the 
relative entropy -F[G|G P ] is nonincreasing, and here we use a quadratic 
approximation of it to show eq. f|23[) . 

We will first prove a needed auxiliary inequality; then, we prove 
eq. d2SD- 



2.1 Some notation 

Here we explain some of our notation conventions, and in particular 
that for weighted LP spaces; for the specific spaces used in our study 
of the linearized coagulation operator, see Section 13.11 

As most of our functions will be defined on (0, oo), we omit the 
limits in integrals when this does not lead to confusion; thus, when we 
omit the limits of integration, it is understood that integration is on 
(0, oo); when we omit the variable, it is understood that integration 
is on the y variable. For example, g(y) dy is understood to be an 
integral over y £ (0, oo), and for a solution g = g(t,y) depending on 
time t and size y, the expression f g denotes J* °° g(t, y) dy. In the same 
way, derivatives are with respect to the size variable unless explicitly 
noted. 

We use several spaces of functions on (0, oo), and we similarly omit 
the explicit mention of the interval. Thus, we use: 

L 1 :=L 1 (0,oo) L\ := L 1 ((0, oo); y k dy) (fc € R). (24) 

In general, for a function y \— > m{y) and 1 < p < oo, we write 

LP(m(y)) := L p ((0, oo); m(y) dy), (25) 

to denote the space of functions / such that y i— > \ f(y)\ p m(y) is in- 
tegrable. Likewise, C and C k denote the spaces of continuous and 
fc-times continuously differentiable functions on (0, oo), respectively; 
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Co, Cg are the corresponding spaces of compactly supported functions 
on (0, oo). 

For any function which is denoted g or h, and which is integrable 
on (e, oo) for any e > 0, we always write G, H to mean the primitives 
given by 

{•CO /'OO 

G(y):= g(x)dx, H{y) :=/ h(x) dx (y > 0). (26) 
Jy Jy 

The same is done for the pair g p , G p , where g p is the stationary solution 
from eq. (|12p . Hence, 

G P (y) := / g P {x) dx = 2e~T (y > 0). (27) 



2.2 Aizenman-Bak's inequality and its linea- 
rization 

The following lemma can be found in [1, Proposition 4.3]. 

Lemma 2.1 (Aizenman-Bak inequality). Take f : (0, oo) — > R with 
f > 0, / G L 1 , and / log / G L 1 . T/ien, 



o jo 



f(x)f(y) log /(x + y) 



< / / f (x)f(y) log f(y)dxdy 
Jo Jo 



roo \ 2 

^ /(s)rfxj , (28) 



and the equality is attained only for the exponential functions f(y) = 
e~w, with fi>0. 

Now we expand this inequality to second order near an exponential 
function to find the following: 

Lemma 2.2 (Linearized Aizenman-Bak inequality). For any function 
h G L 2 ((l + y)e fJ-y ) and any \x > 0, it holds that 

4^ hHeW <2fi(J h\ (J yhj+J^ ^V^+J^J h 2 e' iy ,{29) 
where H was given in \26\). 

Proof. Fix jj, > 0, and take a continuous function h with compact 
support on (0, oo). For small enough e > 0, the function 

f e (y) := e-w + eh(y) (y > 0) (30) 



11 



is in the conditions of Lemma [2, H so Aizenman-Bak's inequality holds 
for the function / e , for all small enough e. One can easily check that 
both sides of the inequality are differentiable in e. As the equality is 
attained at e = 0, it must hold that 



de 2 



r d 2 
e=o de z 



h, (31) 

e=0 



where 1\ and I2 are the left and right hand sides of Aizenman-Bak's 
inequality for f e , respectively. This will give the inequality we are 
interested in, so let us calculate these terms. 

Step 1: First order derivative. Let us calculate the first order 
derivative of all terms in the inequality. The regularity of f e justifies 
the derivation under the integral sign, and we obtain the following for 
small enough e: 

d d f°° f°° 

fe ll = fej J fe( x )fe(y) l °gfe(x + y)dxdy 

f'OO f'OQ 

= 2 1 I h(x)f ( (y) log f e (x + y) dx dy 
Jo Jo 

+ fe(x)f e (y) ffi + j\ dxdy =: I U + J 12 . (32) 

Jo Jo feix + y) 



For the first term in I2, 

1 

fe{x)fe{y) log fe(y) dx dy 



d_ 
~de 



Jo 



"OO poo 

K x ) / fe(y)logf e (y)dydx 
10 Jo 

/'OO /'OO 

+ / fe(x) / h(y)\ogf e {y)dydx 
Jo Jo 

+ ( I h ) ( I fe) =■ hi + I22 + I23. (33) 



And for the last term in I2 we have 

2 



I (f '«<*>*) - s (/*)(/*)-= J - <34> 
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Step 2: Second order derivative. We have: 
hi 



d 



6=0 



2/j, / h(x)h(y)(x + y) dx dy 
Jo Jo 

+ 2 1 I h(x)e^ x h(x + y)dxdy 
Jo Jo 



-4fj, (J ' hj (J yh\ + 2 J hHe m (35) 



d_ 

d7 



poo roo 

I 12 = 2 / h(x)e^ x h(x + y)dxdy 

e=0 Jo JO 

oo />oo 

/ e^ x+y ^h{x + yf dxdy 

JO 

= 2 / /life'* 



y h 2 yeM. (36) 



~de 



6=0 



'21 



e=0 



'22 



-»{J h ) [J " h ) + {J", 



d7 



e=0 



'23 



~dl 



6=0 



'24 



2 



(37) 
(38) 
(39) 
(40) 



Step 3: Final inequality. Putting together eqs. (|35j) and (|36l) 

we have 



de 2 



_h = ^ f hHe m - An (J h^j (J 



yh 



h 2 ye^ y , (41) 



and eqs. ([3"7j) - (fI0|) show that 



de 2 



6=0 



2/i / h 



Finally, writing out the inequality in eq. (|31l) in view of eqs. (|41j) - 
I2j) . we obtain precisely eq. ()29j) . This shows our inequality for any 
h £ Cq(0, oo), and taking a suitable limit of such functions shows the 
inequality for any h £ L 2 ({l + y)e m ) (note that for such an h all terms 
in the inequality are finite). □ 
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2.3 Spectral gap inequality 

Take p > 0, which will be kept fixed in the rest of this section. One 
can linearize the self-similar eq. around the profile g p with mass p 
(given in eq. (I12p ) to obtain dfh = Lh, with L given by the following 
definition: 

Definition 2.3. For a function /i 6 Cq, we set 

Lh(y) := 2h(y) + yh'{y) + 2C(g p , h){y). (43) 

Using the expression of C(g p ,h) from eq. (|10p and noticing that 
h* g p = —^H + g p f h+ - p H * g p (an integration by parts), one has 

2C(g p , h) = h*g p -g p Jh-2h = -2h-~H + ~H* g p , (44) 
so (I43p may be rewritten as 

4 2 

Lh = ydyh - -H + -H * g p . (45) 

Observe that this expression can be generalized to functions h which 
are not necessarily integrable at y = 0, as it is H which appears in 
the convolution; it is for this reason that it will be useful later. The 
following expression for the primitive of the operator C will also be 
needed below: 

Lemma 2.4. For g,h G L 1 , and y > 0, it holds that 

\ ry i r°° 

C(g,h)(x)dx = - / g(x)H(y - x)dx - -H(y) g(x)dx, 
z Jo z Jo 

or, written more compactly, 

C(g,h)=g*H-H j g, 

where H is given in eq. 126\) . 



A proof of the above result can be obtained by direct integration 
of eq. (|10p . and we omit it here. For an explicit proof, see for example 
that of Proposition 9 in |12| . A direct consequence of this is the next 
expression for the primitive of the operator L: 

Lemma 2.5. For h £ Cq and y > 0, the primitive of Lh is given by 

poo 

Ch{y):= Lh(x)dx = -H(y)-yh(y)+g p *H(y), (46) 
Jy 

where H, G p are the primitives of h, g p given in eqs. 1126}) and [UTty , 
respectively. 
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Proof. Using Lemma 12.41 and integrating by parts, 

f'OO f'OO 

Ch(y) = 2H(y) + xti(x)dx + 2 C(h, g p ){x) dx 
Jy Jy 

= H{y) - yh{y) + g p * H -H j f g p = -H(y) - yh(y) + g p *H, 

as the integral of g p is equal to 2. □ 
Proposition 2.6 (Spectral gap inequality). For any h £ Cq such that 

oo 2 poo 

H(y) 2 e~p y dy < oo and / yh(y)dy = 0, 



o 



it holds that 



{h,Lh)_ x ._ <-|WI_i2, (47) 



where (•, ■)_■, a is the scalar product associated to the norm 



-l,2/p> 



r°° 2 
{h,Lh)_ t 2:= / H(y)Ch{y)e— P y dy. 
,p Jo 

Proof. Set n := 2/p as a shorthand. We explicitly calculate this ex- 
pression by using eq. (jlBj) : 

/■OO /'OO 

(h,Lh)_ l4l = J H(y) J Lh{x)dxe m dy 

= - j H 2 e m - J Hhyef*" + J H(g p * H)e m . (48) 
We rewrite these terms using integration by parts: 

J Hhye^ = - j Hhye m + j H 2 e m + fj, J H 2 y , (49) 
and hence 

j I Hhye™ = ~ J H 2 e™ + ~ / H 2 ye m . (50) 
As for the third term in eq. 



poo a poo py A p 

/ H(g p * H)e m = - / H(y) / H(x)e' xx dxdy = - / HHe™ , 
Jo PJo Jo PJ 

(51) 

where Ti is the primitive of H: 

poo 

H(y) := / H(x)dx. 
Jy 
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Putting eqs. (|50|) and (f5Tj) into ([48]) one obtains 

(/i, Lh)_ lfi = ~ J H 2 e^ y - 1 - J H 2 ye> + - p j HHe>. 

Using the linearized Aizenman-Bak inequality (Lemma 12. 2p with H 
instead of h and \x := 2/p gives 

{h ,Lh)_^<-l j -± j H*ye> 

+ - J H 2 yel y + \ J H 2 el y = -J H 2 e> = -\\h\\\^ 

where we have taken into account that J yh = when omitting one 
of the terms in the inequality. □ 



3 Spectral gaps of the linearized oper- 
ator 

3.1 Functional spaces and definition of L 

For \i > and integer k > —1, we consider the following norms for 
a function h : (0, oo) — ► R, defined when the expression below make 
sense: 



2 



poo 

/ y^ k+l \D k h{y)fe^ dy. (52) 
J o 



Here, D k h denotes the k-th derivative of h; for D 1 h we often use the 
more common h', and when k = —1 we set D~ 1 h := —H (cf. eq. (Til 
We also consider the following spaces: 



:= Ih'\H e L 2 {e m ) and J H = j 
X fej/1 := ih | y k+1 D k h e L 2 (e w ) and Jyh = o\ 



(53) 



(54) 



where it is understood that the members of these spaces are at least 
distributions on (0, oo), and the derivatives are meant also in the sense 
of distributions. In each X/.„ for integer k > — 1 and [i > 0, the above 
defined ||-|| fc „ gives a Hilbert norm, and their associated scalar prod- 
ucts are denoted by (-, -) fc . We will prove later that the linearized co- 
agulation operator has a spectral gap in these spaces, and considering 
only functions h with j y h = is tied to the fact that the coagulation 
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equation conserves the mass, so its linearization cannot have a spec- 
tral gap in a space without a restriction of this kind. Observe that 
the restriction J H = in X—i u is the same as J y H' = — J yh = 
when H is regular enough, so it is the same as in the other spaces. 

For given \i > 0, the spaces Xk tfJ , form a scale where each space is 
dense in the next one, as we show in the next two lemmas. 

Lemma 3.1. For any fi>0, any integer n > and h 6 L 2 {y 2n+2 e^ y ) 
it holds that 

/ H 2 y 2n e™ dy<4 /i 2 y 2(n+1) ^ y dy. (55) 
Jo Jo 

Proof. Notice that h is integrable on (e, oo) for any e > 0, so H is 
well-defined and H{y) — > as y — > +oo. We use Cauchy-Schwarz's 
and Hardy's inequality as follows: 

H 2 y 2n e fiy dy = / / (y)h(x)y 2n e m dx dy 

Jo Jy 

POD [X 

= / h(x) / H{y)y 2n e m dydx 
Jo Jo 

/'OO -1 rx 

< \h{x)\x n+1 e% x - \H(y)\y n e2 y 'dydx 



o 



i 

2 \ 2 



< 2 f /°° /*V"+ V x dx^) 2 f /°° H 2 y 2n e m dy 



This proves the inequality in the lemma. □ 

Lemma 3.2. For [i > and integer k > 0, is contained and 

dense in X^i „. 

Proof. Lemma [3.11 proves the inclusion. The claim that they are dense 
inclusions is a consequence of the fact that the set {/ £ Cq°([0,oo)) | 
J yf = 0} is dense in all of them. □ 

Let us also give other useful inequalities: 

Lemma 3.3. For /i > and any h £ L 2 (e m ), 

H 2 e» y < 4 y ^ 2 e ra - (56) 
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Proof. The inequality (2h — fJ.H) 2 > and an integration by parts 
gives 

= 2fi (J h\ + f i 2 J H 2 e™ >fi 2 J H 2 e™. 

□ 

A consequence of Lemma (|56p applied to the function y k+1 D k H 
in the place of H, and Lemma 13. 11 is the following inequality: 

Lemma 3.4. For integer k > and ji > 0, there is a number C = 
C(/c,/x) > such that 



\H\\ k)fl < C \\h\\ Ktl 



for all h € X k 



Another useful property of these spaces is that exchanging deriva- 
tives and powers of y gives other equivalent norms. More precisely: 

Lemma 3.5. Take fi > 0, and nonnegative integers k, a, b, n, m such 
that a + b = n + m = k + \. Assume that h G C k ~ l is such that 

lim y^rrhly) =0 (» = 0, . . . , k - 1). (57) 

Then there is a (constructive) constant K = K(k,a,b,n,m) > such 
that 



y a D k (y h h) 



< K 



y n D k (y m h) 



for any h G L 2 (e^ v ), in the sense that the left hand side is finite 
whenever the right hand side is, and then the inequality holds. 

One can prove this by repeated application of weighted Hardy 
inequalities (cf. [10], Example 0.3]), and we omit the details. 

Remark 3.6. Prom this it is easy to see that the L 2 (e ra )-norm of any 
combination of products by positive powers of y and derivatives of h 
gives an equivalent norm as long as the total sum of all the powers 
and that of the orders of all derivatives is fixed. It is important to 
restrict these considerations to functions h satisfying (|57p and having 
some integrability property at +oo, (e.g., h G L 2 {e m ) in our case), as 
otherwise the finiteness of the right hand term does not imply that of 
the left one. 
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The above equivalence of norms can be used in our X^^ spaces for 
H, k > 0, as the condition in the previous lemma holds: 

Lemma 3.7. If h £ Xj.„ with fi > and integer k > 0, then 

lim y t+2 D l h(y) = (i = 0, . . . , k - 1). (58) 

Proof. Take an integer < i < k. As h is in C one 

directly sees that the first derivative of y' l+2 D' l h(y) is integrable near 
y = 0. Hence, y t+2 D t h(y) has a limit at y = 0. As y l+1 D l h{y) is 
integrable, it must be that this limit is 0. □ 

One can define the operator L in the spaces X^^, of course agreeing 
with our previous one in Definition 12,31 whenever both are applicable. 
For the rest of Section [3] we fix p > 0, and consider the linearization 
of the nonlinear operator in eq. © around the self-similar profile g p : 

Definition 3.8. Take /i > 0. For h £ Xq^ (in particular, for h £ Xj. ^ 
with k > 0) we define Lh by eq. (|45j) . where the derivative of /i is 
taken in the sense of distributions. 

Our aim in this section is to show that the operator L is defined 
as L : Xk+i n — ► -Xjfc,« for k > —1, and that it is a closed operator in 
the spaces X^^ for fx > and any integer fe > — 1. Before proving 
this we will need two lemmas in which we study the operator C{g, h) 
which appears in the definition of L. 

Lemma 3.9. For ji > 0, integer k > 1 and g,h S X^ it holds that 

D k {y k+ \g *h)) = Y,( k + X ) (D k+1 - l (y k+1 - l g)) * (D^tfh)). 
i=o ^ ^ ' 

In /act, t/iis ZioWs w/ien y k D k+1 g,y k+l D k h G L 2 (e ra ); i/iis is, w/ien 
<7, /i satisfy the conditions for being in X^n, but without imposing that 
J yh = J yg = 0. In particular, this applies to g = g p and h £ X^ 4l . 

Proof. One has 

rv 

y k+1 (g * h){y) = / y k+l g(x)h(y - x) dx, 
Jo 

and the identity in the lemma can be obtained by writing the binomial 
expansion for y k+1 = (x + (y — x)) k+1 and differentiating the resulting 
convolutions. The border terms which appear when differentiating the 
integrals between and y vanish due to (|58l) . □ 
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Lemma 3.10. ForO < fj, < 4/p and integer k > —1 i/iere is a constant 
K = K(p, k, p) such that 

\\C(h,g p )\\ k ^<K\\h\\ k>ll forallheX k>tM . 

Proof. Here, C(h,g p ) is understood to be defined by eq. (PHI) , which 
makes sense for /i £ ^A;,^ (which implies if G L 2 (e^ y )). Numbers 
K\,K2, ■ ■ ■ below are assumed to depend only on p, k and p. 
Let us first prove the lemma for k > 0. From 

2C(g p ,h) = -2h-~H + ~H*g p 

and our previous lemma, we deduce that 

2 \\C(g p , h)\\ ktfM < 2 \\h\\ ktll + ^ + - p \\H * ^|| fc>|t . (59) 



The second term in eq. (i59j) is bounded by K\ 
31 for some K\ > 0. For the third term in (|59|) we have: 



thanks to Lemma 



\H*g p \\ <K 2 \\D k (y k+1 (H*g p ) 



k+l 



< 



i=0 



L 2 ( e w) 

D^Ky^gp) 



< K 4 \\h\ 



for some constants K 2 , K3, K4 > 0, where we have used Lemma 1331 in 
the first inequality (observe that the conditions at are met), Lemma 
3.91 in the second, and Lemma 13. II for the third one. As the integrals 
bounded in the last inequality are finite as long as < p < 4/p, the 
lemma is proved for k > 0. 

Now, for k = — 1, by density it is enough to prove the inequality 
when h is Cq°([0, 00)). We use the expression of the primitive of C(g, h) 
from Lemma 12.41 to obtain, using Young's inequality, that 



9p 



nC{g p ,h)\\-x^< \\9p*H\\me^y) + \\ h \\-i, P J: 
< \\H\\l^) j \9 P \ei y + \\h\\_^j g P <2\\h\\_ l4i j \g p \e^. (60) 



□ 



We finally have the following: 
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Proposition 3.11. For any integer k > — 1 and < fx < 4/p, the 
operator L is defined between X^+i^ and -Xfe,^-' 

L ■ Xk+\^ — > -^fe,/i- 

5een as an unbounded operator on X^^ for integer k > — 1 ana 1 reaZ 
fi > 0, L is a closed operator with dense domain. In addition, for 
every [X > and integer k > — 1 iftere is a constant K = K(p,k,fj,) 
such that 

(Lh,h) k ^<K\\h\\ 2 k ^ (/!£%,,), (61) 

and consequently L generates an evolution semigroup in each of the 
spaces Xk 4l . 

Proof. To see that L is denned between X^+i^ and „ for any k > 
— 1, take h £ X^+i^ and notice that every term in (|%3"1) is in Xfc jAt , as 
can be seen from Lemma 13. 101 and the fact that h £ -Xjfc+i,» implies 
that yh! £ X^^ (see Lemma 13.51 and the remark that follows). 

For k > —1, let us prove that L : X^ + i^ — > is a closed 

operator in the space Xfc^. Take a sequence ft n in Xk+i^ which 
converges to some h £ X/~n in the norm of „, and such that Lh n — > 
ft. in Xk tfl for some ft £ X^^. We need to show that ft G Xfe + i 5 ^ and 
Lh = ft. For such a sequence we have 

Lh n = 2h n + yh' n + 2C(g p , h n ) 

As the sequence {ft n } converges to ft, it is clear from Lemma 13.101 
that the first and last terms in the expression of Lh n converge to 2ft 
and 2C(g p ,h), respectively, in X k ^. As Lh n converges to ft, it follows 
that the second term, yh' n , converges to something in X^^. From 
the equivalence of norms in Lemma 13.51 this implies that in fact h n 
converges to something in the norm of Xj-+\ ^. As this is a stronger 
norm than ||-|| fe , we see that the limit must be ft, which implies that 
ft £ Xk+i^ and Lh n — > Lh in X k n. Hence, ft = Lh. 

Finally, let us prove the inequality (|6ip for k > (for k = — 1 a 
similar argument proves it). From expression (|43|) we have 

(Lh, h) k>fi = 2 \\h\\l tfM + (yh', ft)^ + 2 (C(g p , ft), ft)^ 

< 2 \\h\\%, + W> h )k,» + 2 ^)H^ H%,„ • (62) 

For the second term, as D k (yh') = kD k h + yD k+1 h, 

(yh',h) k ^ = k\\h\\ 2 k ^ + J y 2k+3 (D k+1 h)(D k h)e»y < k\\h\\ 2 kfi , (63) 
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as an integration by parts shows that the term we omitted is negative: 
J y 2k+3 (D k+1 h){D k h)e^ y 

y 2k+2 {D k h) 2 e^ - | J y 2k+3 (D k h) 2 e^ 

= -^\\h\\l»-^Jy 2k+ \D k hfe»y. (64) 

For the last term in (I62p . using Lemma 13, 101 

\\C(g P ,h)\\ k ^<K\\h\\^ (65) 
Hence, we finally obtain 

<LM> fciM <iriiifciiiU 

for some K\ depending on p and p. □ 

3.2 Spectral gap in X-\^i p 

A direct consequence of Proposition 12.61 and a limit argument is the 
exponential decay of the evolution semigroup defined by L: 

Lemma 3.12. For h G X 02 / P , 

(h,Lh)_ 1>2/p <-\\h\\ 2 _ 1>2/p . 
As a consequence, for h G •X'— 1 2/p we have 

ll e ^ ||- W ^ll fc0 Li >2 /p e " t (*^°)- ( 66 ) 

3.3 Extension of spectral gaps to lower expo- 
nential weights 

The evolution semigroup generated by L in the spaces Xf. „ has the 
remarkable property of creating exponential moments of (D k h) 2 in 
finite time. This will allow us to prove that, if L has a spectral gap in 
X/i V for some k > — 1, then it also has a spectral gap in Xk tP for any 
< p < v. 

Lemma 3.13. Take < p < v < 4/p and an integer k > —1. Ifh is a 

solution of the linear self-similar equation dth = Lh with h(Q) £ X^n, 
then the norm \\h(t)\\ k v is finite at time t = to := \og(u/ p), and 

l|fcft>)IU,„< Q \W)\\k* 

for some positive constant K = K(p, k, p, v). 
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Proof. Let us first prove the lemma for k > 0. From eqs. (|62p - (|65p in 
the proof of Proposition 13.111 for any < 7 < 4/p we have 

(h, Lh) Kl < K \\h\\l^ -lj y 2k+3 (D k h) 2 e™, (67) 

where the constant K = K(p, k, 7) > is the same as that of Propo- 
sition 13.111 With this, take (fi(t) := pe t and carry out the following 
computation: 

^\\h\\l m = 2(h,Lh) km +<p'(t) [ y 2k ^(D k h) 2 e^ 



alt 

< 2K(d>(t)) \\hf m + (<P>(t) - m J y 2k + 3 (D k h) 2 e™ 

= 2K(<t>(t))\\hf m , (68) 

where K(<j>(t)) is the constant obtained from (I67p . for which we write 
the dependence on </>(i) explicitly. Now, for t < to, 

I \\h\\l it) <2K(m\\h\\l {t) < Ki\\h\\ 2 m (69) 

for some K\ > 0, as the constant K((f>(t)) is bounded for < t < to, 
which can be checked from the proof of Proposition 13.111 Hence, as 
0(0) = M, 

.2 ^ Jdtonou* - f"\ Kl \\hO\\ 2 



which proves the lemma for k > 0, given that (f>(to) = pe to = v. 

For k = — 1 the same argument can be carried out by using the 
inequality 

(h, Lh) < K \\h\\\ y - 1 f yH 2 e™ 
instead of (fBTl) , and the following instead of (f68l) : 

j t \\h\t mt) = 2(h,Lh)_ im+( p'(t) J ye^H 2 . 

□ 

Lemma 3.14. Take < v < 4/p and an integer k > —1. Assume 
that the operator L has a spectral gap in X^ u of size 5 > 0; this is, 
there exists C > 1 suc/i i/iaf 

lle^lL <C\\h°\l e~ st . 

II II K,J/ — II II K,^ 

Then, L has a spectral gap of the same size in « /or any < p < v; 
this is, there is C = C'(p, C) > 1 such that 

\\e tL h°\\ h <C'\\h°\\ h e- 5t . 

II II k,fi — II II k,fi 
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Proof. Take h in the domain of L as an operator on X^^ (i.e., X^+i 
if > 0, and L 2 (ye m ) if k = —1). Consider the solution t i— ► h(t) 
to eq. ([8]) with initial condition h°. The idea is to estimate ||/i(*)llfc ^ 
using eq. (jBTI) until a time to for which ||/t(£o)||_i v is finite, and after 
that time use the fact that the latter norm is exponentially decreasing 
thanks to the spectral gap we are assuming in the norm of v . 
Choose to := log jj- Then, for t < to, the estimate in eq. (joTj) gives 

j t \\h(t)\\l^ = 2(h(t),Lh(t)}^<K 1 \\h(t)\\l^. 
So, for < t < to we have 



. K 

V 



At t = to we can use Lemma 13.131 with v := 2/p to get that 

IIM*o)|ll,, < H^ll^ K/i 

where K is the one from the lemma. Now, as we are assuming a 
spectral gap of a certain size 5 in Xj cu , we have for t > to 

i,»<\\hml„<\Wo)f k>u e- 2s ^ 



J 1 

Together with eq. (|7(J|) . this shows that 



K / \ K+28 



-) 



<\\h"\\: (-) e - 2a (*-*o) = \\h°\\1 {-) e~ 2St 



— II \\k,n 

with 

C 2 : = max 



<C\\h°\\ k e~ 5t forallt>0, (71) 



, K+28 / \ K x 



The result is extended to all h° G X^^ by density. □ 

Then, from the spectral gap we showed in Lemma 13.121 we obtain 
the following: 

Corollary 3.15. Take < p < 2/p, and consider the semigroup e tL 
generated by L in the space X-\ ifl . Then there is a number K = 
K{p) > 1 such that 

||e' L /i°||_ 1>M <K||/i°||_ liAt e-< for any h° G X_^. 
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3.4 Extension of the spectral gap to Xk :fl with 



Let us shortly discuss the strategy of the proofs in this section. One 
could try the following argument in order to obtain a spectral gap for 
a norm involving a certain derivative D k h of order k > 0. We could 
hope that, when calculating 



one obtains a negative multiple of this same norm, plus some terms 
which involve L 2 norms of h with some weight, and which can already 
be controlled by the semigroup decay results in previous sections; this 
idea can be arrived at by trying it for k = 0, in which case it works 
perfectly. But for k > 0, the term which involves the same norm is 
actually a positive multiple of it, so our tentative argument does not 
directly give a spectral gap in this case. It is for this reason that we 
consider instead norms with a power weight y m for a certain m: one 
is forced to raise the power m if one wants to obtain a spectral gap 
involving L 2 norms of higher derivatives of h. Then, another problem 
appears: if one calculates, say, the time derivative of J y m (h') 2 e liy , 
then some terms which involve L 2 norms of h will appear with the same 
weight y m as we were using in the first place: this is, a term involving 
J h 2 y m e^ y will appear, for which we have no previous spectral gap 
results. Hence, to close the estimates we use a spectral gap in a norm 
with an exponential weight of order slightly higher than the one we 
started with. 

Lemma 3.16. Take < fx < 2/p and an integer k > 0, and consider 
the operator L in the space X^n, with domain X^+i^- Then, L has 
a spectral gap in X^^; more precisely, there is some constant K = 
K(k,fJ,,p) > such that 



for any h° £ X kifi . 

Proof. Take k > 0, assume the result true up to k — 1 (the case k = — 1 
was proved in Corollary 13. 15[) . Denote by h the solution of the linear 
equation with initial condition h° (i.e., h(t, •) = e tL h°), and in order 
to simplify the notation write ||-|| instead of ||-||j.„ and (•,•) instead 
of (•, -) k . The numbers K, K\,K<i, ■ ■ ■ which appear in this proof are 
understood to be positive and depend only on k, fj,, p and 5. 

Pick fj, < v < 2/p. Then, after the time io given in Lemma 13.131 
the norm „ is finite and we have, from the expression for L in 



k > 
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eq. (05]): 

ld_ 

2dt 1 



2 = (h, Lh) = (h, yh') -t(h,H) + - (h, H * g p ) 



P P 

For the first term in ([73]) . eqs. (f6"3j) - (f64"|) show that 



[h,yh!) < -- 



(73) 



(74) 



For the middle term in eq. ([73 



(h,H) 



{D k h){D k H)y 2k+2 e m dy 

(D k h)(D k - 1 h)y 2k+2 e^dy 
--(k + 1) J (D^hfy^e™ dy 

+ ^ J \ D k-l h) 2 y 2k+2^y dy 



fe— l.i/ ' 



(75) 



where in the case k = we denote D~ 1 h := — iT. Observe that 
here the exponential weight has been changed from \i to v, as the 
power weights in the norms which appear are not the right ones for 



the 



I As— 



norm, but higher. We overcome this difficulty by using a 



norm with a slightly higher exponential weight. 

To bound the middle term in (I73p . we will need the following ex- 
pression for D k (h * g p ), which can easily be proved by induction: 



fc-i 



i=0 



D k (H *g p ) = H* (D k g p ) + ^(D^h^D^g^O). (76) 

From this, 

\\H*g p \\= y k+1 D k (H*g p )\ 



< 



k+1 (H*(D k g p )) 



< 



L 2 (ew) 
k—1 



L 2 ( e ™) 

y k +\H*{D k g p )) 



y k D i - 1 h 



+ K 2 \\h\ 



L 2 (eW) 

fc-1,,, (77) 



where the last bound is possible because all the norms inside the sum 
are bounded by a constant times ||^|| fc using a higher exponential 
weight as before. For the remaining term in (|77|) . 

fe+i 

y k +\H * (D k g p )) = £ (* + 1 ) tfB) * (y k+1 ' l D k g p )) 

i=0 ^ 1 ' 
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and hence 



y k+1 (H*(D k g p )) 
fe+i 



< 



< 



E 

i=0 
fe+1 

E 

i=0 



L 2 ( e ™) 
fc + 1 



(y'tf) * (y k+1 - l D k g p )) 



fc + 1 



<^3 



fe+l-i 



L 2 (eW) 



D k g P 



-l,u 



< Ka 



k-l,v ■ 



Notice that the integrals involving g p are numbers that depend only 
on k, fJL and p. Prom ([77]) . 



|# *<?p|| <^5 



lfc-l,f ' 



so 



|</i, if * g p )\ < \\h\\ \\H * g p \\ < K 5 \\h\\ \\h\\ k _ lyU . (78) 

Finally, all three terms from (|73[) are bounded in (|74f) . (I75p and (|78p . 
and together give 



1 d 
2dt 



2 < -| H/if + -i^i H/illLi,, + Z K 5 \\h\\ \\h\\ k _. 



P 



2 
P 

<-^\\hf + K ( 



6 \\ n \\k-lM > 



where we used the elementary inequality ||/t|| fe _ 1 w < e 2 ||/i|| 2 + 
(1/e) ||fo||f_i„ for an arbitrary e > 0. Thanks to the spectral gap in 
IHIfc-i v anc ^ Lemma 13.131 we know that for any 5 < 1 there is some 
constant K > such that \\h(t)\\ k _ hu < K \\h(t )\\ k _ h „ e _<5 (* _ *°) for 
t > to. Hence, for t >to, 



k-l,u 



-28(t-t ) 

< -\\\hf + K^f 



A Gronwall lemma, together with inequality (|6ip for t < to, finishes 
the proof. □ 



4 Behavior of moments 

For later use we will need to have some information on the time evolu- 
tion of exponential moments of solutions to the coagulation equation 
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([8]), which is also interesting by itself, as this equation is one of the 
few particular cases where one can find an explicit expression for the 
evolution of moments. For a function g : (0, oo) — > R we denote 

/•oo 

M k [g] := / y k \g(y)\ dy for k G R, (79) 
J o 

/"OO 

EM := / e^| 5 (y)| for M G R. (80) 
■/ o 

In general it is slightly simpler to solve explicitly the evolution of 
moments in eq. (pTJ) than in eq. ©; but of course, one of them being 
a rescaling of the other, there is a simple relationship between the 
moments of their solutions: if / and g are related by the change of 
variables in eqs. (J6|)-(J7|), then the relationship between moments of / 
and moments of g is 

M k [g(t)]=e t ^M k [f(e t -l)} (81) 
M k [f{t)} = {t + l) k - l M k [g{\og(l + t))}. (82) 

That between exponential moments of / and g is 

E fl [g(t)) = e t E lie - t [f[e t -l)}, (83) 

E^[f(t)] = T ^ [ E IM{t+1) [g(log(t + l))}. (84) 

The evolution of moments of order 0, 1 and 2 for a nonnegative 
solution / to eq. ([I]) (or the corresponding solution g to eq. (JSj) ) is 
easily obtained and well-known: the first moment is a constant, called 
its mass. For the moment of order 0, 



so 



MM = t+ 2/M 0[ /°r <85> 

Using the relation in eq. (I8ip (with A = 0) we have 

M *l = l- e - + 2 2 e-/M„b T < 86 > 
The second moment of / in eq. (H|) with a = 1 satisfies the equation 

where p is the mass of the solution. Hence, 

M 2 [f]=M 2 [f] + \p\ (87) 
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and the second moment for the self-similar equation is 

M 2 [g] = e-'M 2 [ 5 °] + ^p 2 (l - e"*). (88) 

For exponential moments one can also find an explicit expression: 
again from eq. ([I]), and using eq. (|85|) . 

±E,[f}= 1 -E,[ff-M [f}E,[f] 
= 2 E ^\f\ 2 ~ t + K E ^[f], 
with K := 2/Mo[/ ]. This has an explicit solution: 

W*)] = ttV+ 2 2 -I . (89) 

where M ° and £° denote M Q [/(0)] and ^[/(O)], resp. 

Using eq. (1830 we obtain the evolution of exponential moments for 
the self-similar equation with A = 0: 

EpbQogis))] = ^-^ + g 2 5 — • (90) 

This directly implies the following lemma: 

Lemma 4.1. If g is a solution of eq. (|8|) with initial data g° and 
fi > 0, then En\g\ is finite for all t > if and only if the initial data 
satisfies 



2 

1 



E° e - M ° < -n for all0<9<n. (91) 



E^[g] is uniformly bounded for all t > if and only if the initial data 
satisfies the following for some v > /i: 

E° e - M ° < for all0<9 <fi. (92) 

In particular, if E u [g] is finite for all t > 0, then for each /i < v the 
moment E^[g] is uniformly bounded for all t > 0. 

Next we prove that if some positive exponential moment is initially 



2 
P 

after some time. The precise result is the following: 



finite, then every exponential moment less than - becomes bounded 
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Lemma 4.2. Let g be a solution of eq. with initial data g° . As- 
sume that there exists \x > such that E^g ] < oo. Then, for each 
< v < | there is a time T v > (which depends on the initial condi- 
tion g°) such that 



E u [g(t)] < oo for all t > T u . 



(93) 



As a consequence, for each < v < - there is a time T* > and a 
constant K u > such that 

E v [g(t))<K v forallt>T* v . (94) 



Proof. The uniform bound eq. (|M|) is a consequence of eq. (|93|) in 
view of Lemma 14.1} so we will only prove eq. (I93p . Take < f < ^ . 
In view of eq. (|90|) . it is enough to show that, for some /j, > e > 0, 

2 



^b ]-M b°]< 



S - 1 



for all 6 < e, 



(95) 



and then eq. (|93l) holds with T v := log^. As eq. (1951) is a local 
inequality, it is enough to study both terms near 9 = 0: they are at 
6 = 0, and their derivatives are: 



^9 I > !{!n<H,i,h) 



d 



0=0 



yg{y) dy = p 



Hence, the inequality 



2v 



0=0 



0=0 



2 
v 



2 

v < - 
P 



shows that eq. (|95l) holds for some e > 0, and the result is proved. □ 



5 Local exponential convergence 

In this section we will prove Theorem [T31 which is a local exponential 
convergence result of the following kind: we show that if a solution g 
to eq. ([8]) is initially close enough to the stationary solution, then it 
converges exponentially fast to it with a given rate in the norm ||-|| fc 
for < n < 2/p and integer k > —1, as long as this norm is finite at 
time t = 0. 

Here, the concept of close enough is measured in the sense that the 
initial relative entropy to the equilibrium is small. The reason for this 
is that this closeness must be enough to guarantee that the nonlinear 
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term in the evolution equation is small when compared to the linear 
part, which we know is well-behaved. However, if g is a solution, the 
smallness of the nonlinear part essentially involves norms of the kind 
f \g — g P \, which cannot be controlled by our norm \\g — g p \\_ 1 2 / p - 
Consequently, we need to measure closeness in some way that is prop- 
agated in time: what is needed in the proof of the theorem is that, for 
e > small enough, 

poo 

/ \g{t,v) - g P {y)\ dy < e for all t > 0. (96) 
J o 

We do not know of any simple condition on the initial data that guar- 
antees this to hold except requiring the relative entropy to the equi- 
librium g p to be small. Then, as the L 1 norm of g — g p is controlled by 
the relative entropy, we know that (|96[) also holds. When we consider 
convergence in norms ||-|| fc „ with k > we could control f\g — g p \ 
by the norm ||<7 — g p \\ k and obtain a result where closeness to g p is 
measured in terms of \\g — g p \\ k „, but this is less convenient for the 
global result in Section El Because of this, we always state our results 
with an entropy condition. 

Before proving our local convergence theorems we will need two 
previous lemmas. In the first one we show that the function _F[g|<7p] 
controls the distance of g to the equilibrium g p in the L 1 norm, and the 
second one will be used to estimate the nonlinear term C(g — g p , g— g p ) 
which will appear when calculating the time evolution of \\g — g P \\ k ■ 

Lemma 5.1. There is some constant K > such that 



coo 



J \<lUl) - <///< A" max F[g\g P ]. F >|ry,,] 



Proof. We can write 



where ^(y) := (y + 1) log(y + 1) — y > 0. As it happens that ^ is a 
convex function for which \&(y) > ^(|y|), from Jensen's inequality we 
have 



9(y) l 



9 P (y) 



dy 



poo 

F[g\g p ) > / g p (y)V 
Jo 

>2*(±J \g{y) - g p {y)\ dy\ 
where the 2 appears from g p = 2. Hence, 

J \g-g p \<2*- 1 (^F[g\g p ] \ < if max L] F[g\g p \, F[g\g p ]\ 
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for some constant K > 0, which proves the lemma. 



□ 



Lemma 5.2. For p > and < fx < A/p there is a constant K 
K(p) > such that 



\\C(h,h)\\_^<K\\h\\_^l \h\e&dy (hex. 



(97) 



Also, for integer k > and p < v < A/p , there is a constant K 
K(k,p,v) > such that 



\\C(h,h)\\k*<K\\h\\ k>lt J \h\e&dy 



+ K 



k-l,v 



(98) 



for all h G Xk v . 



Proof. Eq. (|9"7]) is a consequence of the calculation in eq. (foT)]) with 
instead of g p . Now, for k > 1 and /i G -Xfc ; „, from the expression 
C(h,h) = ^h * h — h J h we have 



l|tf(M)[|^< 5 l|fc*%^%,, 



(99) 



and from Lemmas 13.51 and 13.9 



\h * 



fc+1 

i=0 
< ^2 



L 2 ( e ™) 

(D i ^ 1 (y i h)) * (D^^^y^ 1 -^)) 



L 2 ( e ™) 



i=2 



fe,/U 



k.fi 

k + 1 



J \h\e2 y + K : 



5 l^llfc,^ / \h\ye 

D k+l-i^ y k+l-i h ^ 



j \h\e* v + f \h\yezA + K, 



eNLi,,- (100) 



Here, we have bounded the terms in the sum using Young's inequality; 
the terms for i = 0, i = 1 and i = k + 1 have been bounded putting 
the 1? norm on the part with k derivatives (and an integration by 
parts for i = 0). The rest of the terms have been also bounded using 
Young's inequality, and we have applied Cauchy-Schwarz for the last 



inequality, together with the fact that \\h\ 



k-l,u 



controls all norms 



v with < i < k. Notice that K\ through depend only on k 
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and fj,, and Kq only on and v. Finally, from (|99p . 



\\c{Kh)\\^<\H k , 

<K 7 \\h\\ k ^f \h\e& + K t 



6lHlfc-l,* 



2 

fc-l,i/ • 



This proves our result for k > 1. For the case = 0, the same 
calculation proves the inequality, with the difference that in this case 
the term multiplying K4 in (|1UU|) does not appear and the argument 
is much simpler. □ 

We can finally prove Theorem ll.4t 



Proof of Thorem \1.4\ Step 1: Proof for k = —1. Denote ||-||— 1 ^ 
by ||-||. We assume that \\g° — g p \\ < 00, and want to show the expo- 
nential convergence in the same norm. We prove an a priori estimate 
below, which can be made rigorous by considering approximate reg- 
ular solutions instead of solutions to the full equation. The numbers 
K,Ki,K2, ■ ■ ■ which appear in the proof are understood to depend 
only on p, E, 5, k, \x and v. 

Consider a solution g of the self-similar equation eq. (JSj) in the 
conditions of the theorem, and set h := g — g p , h° := g° — g p . Then h 
is a solution of the equation dth = L(h) + C(h, h), and by Duhamel's 
formula we have 



h{t)=e tL h°+ e^ L [C(h(s),h{s))]ds (t > 0). (101) 



Using our bound on the exponential decay of the semigroup e tL from 
Corollary 13 . 1 5 1 we obtain 

< K x e~* + K x f \\C(h(s),h(s))\\ e-( t ~ s '> ds, 
Jo 

for some K\ > which depends on p, E and /x. Now, use the bound 
in Lemma 15.21 to get 

<K l \\h Q \\e- t + K 2 I EE[h{s)}\\h{s)\\e- {t ~ s) ds, (102) 
Jo 2 

where En[h(s)] is the exponential moment of order p,/2 of h (defined 
in eq. (|80|) ). Considering E and v from the hypotheses of the theorem, 
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we have E u / 2 [g(t)] — E f° r an t >0, so there is some constant C\ > 
such that 

E u/2 [h{t))<Ci for alii >0, 
(as h = g — g p , and g p {y)e^ v is integrable). By Holder's inequality, 

E±[h(s)\ < M [h(s)f^ E u/2 [h(s)]v < e 1 C 2 =: e 2 for s > 0, (103) 
where 

C 2 :=C{, ef 71 := Kmax{V~e,e} , 

and K is the one from Lemma 15. 11 We have used here the bound on 
the relative entropy in the hypotheses, and the fact that it is conserved 
in time, as the relative entropy is nonincr easing. 
Using this bound on eq. (|102p we have 

\\h\\ < K x \hH e-* + e 3 f \\h{s)\\ e-(*- s ) ds, (104) 

Jo 

where e 3 := €2-^2- Equivalently, 

INe^tfill^H+es f \\h{s)\\ e s ds, 
Jo 

and by Gronwall's lemma, 

||/i[|e* < K x ||/i°||e £3t , 

or 

\\h\\ < K^W^We-t 1 -^. 

Choosing e so that 1 — 63 > 5 proves the theorem. 

Step 2: Proof for k > 0. For this step, denote by ||-|| the 
norm . The difference in this case is that the inequality (I98p has 
an additional term, which can be dealt with by using our previous 
steps. Take k > 0, and assume that we have proved the result on 
the convergence in || - ||^,_^_ for < y < 2/p. Repeating the same 
reasoning as before we arrive at the following instead of (I104p : 

\\h(t)\\ < K, \\h°\\ e-* + e 3 f \\h(s)\\ ds 

Jo 

+ K 2 f\\h{s)\\l_^e-^ds, (105) 
Jo 

where the inequality (|98p was applied with some v 1 such that /j, < v' < 
v in the place of v. Now, as we know that ||/i(s)|| fc _ 1 is exponentially 
decaying in time, for any 5 < 5' < 1 we have 

\\h(a)\\ h _^ <K 3 ||^|| fc _ v e- 5s . 
34 



Hence, we get 

||h(t)|| < Ki ||/i°||e-' 

+ 6 3 f \\h{s)\\ ds + K A v e' 5 ' 1 

Jo 

< K 5 \\h°\\ k e- 5 ' 1 + 6 3 f \\h{s)\\ e-*'(*-) ds, 

Jo 

using that controls both H^H^-i^ and Now this in- 

equality has the same form as (I104|) . and the same reasoning followed 
in the first step finishes the proof. □ 



6 Global convergence 

To prove global convergence we consider two stages in the evolution of 
the equation: when the solution is close enough to equilibrium, we can 
use Theorem 1 1.41 to get an exponential rate of convergence; when it is 
far from equilibrium we use the following lemma, which says that the 
solution converges to equilibrium exponentially fast in the L 2 norm: 

Lemma 6.1. Take a nonnegative initial condition g$ £ L\, with 
derivative g' G L 1 and mass p > 0, and consider the solution g to 
eq. (|SJ). There are explicit constants C, T > which depend only on 
go such that 

\\g{t)-g p \\ 2 <Ce-¥ (t>T). (106) 

The lemma is independent from the rest of the paper, as it can 
be proved by using the explicit solution of the Fourier transform of 
equation (jSJ). The technique to prove it is very similar to the one used 
in [TS] to obtain the uniform convergence to equilibrium; the main 
difference is that here we are interested in obtaining a rate, while 
in [15] emphasis was placed in giving very general conditions on the 
initial data under which convergence holds, without any given rate. 

Lemma 16.11 does not give the optimal rate of convergence, but 
we only need it to show that any solution, after some initial time 
which may be explicitly given in terms of the initial datum, is in the 
conditions of our local result in Theorem 11.41 an d hence prove our 
global result in Theorem 11.61 However, one could be more careful in 
the proof below and, imposing slightly stronger decay conditions on 
the initial datum go> obtain the optimal rate e~ l . 

We will first prove Theorem 11.61 assuming Lemma 16.11 leaving the 
proof of the lemma itself for the end. 
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Proof of Theorem \l.b\ Let us show that, for any p < u, we can apply 
Theorem 11.41 for the exponential decay of the norm after an 

certain initial time to- 

Lemma 14.11 and an argument as in the proof of Lemma 14.21 show 
that for some < p! < u, the exponential moment £y /2 [g] is uni- 
formly bounded for all times t > 0. So, take < p < p! with p < 2/ p. 
Then, the first condition in Theorem 11.41 holds. Together with this, 
Lemma 16.11 shows that the relative entropy F^lg^] tends to zero at 
an explicit rate, and hence condition 2 in Theorem ll.4l is also satisfied 
after some initial time to- 

Finally, the same calculation as in Theorem 11.41 shows that, for 
< p < p! with p < 2/p, the norm ||<7 — gpW^. is bounded in bounded 
time intervals, and in particular it is finite at t = to- This is enough 
to apply the Theorem 11.4^ so we have 

Wg^-gpW^^K^ 1 -^ (t>t ), 

for any < p, < v with p < 2/p. As \\g(t) — <?p|| fc „ is bounded for 
t S [0, to], this proves the result. □ 

Proof of Lemma \6.1\ By the scaling and time-translation symmetries 
of eq. ([Tj) , it is enough to prove the result when J y g = J g = 2 (see, 
e.g., [IS]). We consider the Fourier transform of g, 

«MM):= / e-*Wg(t,v)dy. (107) 
Jo 

A very similar calculation to that in section H] shows that (f> has the 
following explicit expression, identical to (|90p (notice that we are as- 
suming / g = Mq = 2): 

_|__ T + 1 2+(t-1)(2-0o(£)) 

for p E M. and t > 0, denoting r := e*. By Plancherel's theorem, 
proving our result is equivalent to showing that 



J oo || 2 



< Ce"5' (t > 0), (109) 



where c^oo is the Fourier transform of the equilibrium with mass 2, 
g 2 (y) := 2e~\ 

<t>M ■= r-^r- (ji G R). (HO) 
1 + Z/i 
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Step 1: Approximations of <fio(x) at x = As 0(0) = 2 and 
d 



dx 



(j) (x) 



d_ 

dx 



(j>o(x) 



oo 

-i / e- ixy yg (y)dy, 
lo 

/•oo 

< / ygo(y)dy = 2, 

Jo 



we readily have 

\4>o(x) ~ 2| < 2|x| E 
Similarly, 0'(O) = — 2z and 
d 2 



<ix 2 



0o (a;) 



ia V5o(2/)%, 



dx^ o( ^ 



< / y 2 go{y) dy =: M 2 , 
'o 



and then, using Taylor's series at x = 0, 



x) - 2 



+ 2i 



< M 



(x E 



(111) 
(112) 

(113) 

(114) 
(115) 

(116) 



Step 2: Lower bound of the denominator in ( 11081) . For 

and r > 1 we have, calling x := h/t, 

|2 + (r - 1)(2 - o (x))| > |2 + 2* M | - |2»> - (r - 1)(2 - ^o(x))| , 

(117) 

and rewriting this second term gives 

2i/i - (r - 1)(2 - o (x)) = M ( 2< + ^° (ic) " 2 > ) + 2 - o (x), 



/z/t 



and using eqs. ()113f) and fjl 16j) . 



\2ifx - (r - 1)(2 - o (x))| < ^M 2 |x| + 2 |x| (x = — E 

Z T 

and so from (11171) . 



(118) 



|2 + (r - 1)(2 - o (x))| > 1 + H - -M 2 \x\ - 2 |x| . 
Finally, for |x| < ei := min{l/4, 1/M2}, 

|2 + (r - 1)(2 - o (x))| > 1(1 + | M |) (|x| = H < ei ). (119) 

z r 
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When \x\ > e±, we estimate the denominator the other way round. 
Call 

e 2 := inf \2 - <fo(x)\ > 0. 

|x|>ei 

The above defined e 2 is strictly positive because (f>o is continuous, 
|(/>o(x)| < 2 for \x\ > and 



|0o(^)| < 2d— (|x|>0), 



(120) 



with 



Ci := 



Wo(y)\ dy. 



Hence, for \x\ > e\ the denominator can also be estimated as 

|2 + (r - 1)(2 - o (x))| > ( r - 1) |(2 - fa{x))\ - 2 > (r - l)e 2 - 2, 

(121) 

and then for r — 1 > C 2 := 3/e 2 , 



|2 + (r-l)(2-^o(x))| > 1 (|x| 



Putting together (jll9[) and A122 jt we obtain 



>ei, r>C 2 + l). 

(122) 



|2 + (r-l)(2-^,( M /r))|>- (/x G R, r > C 2 + 1). (123) 

Step 3: Estimate for large. Take K > 0. We have, for 
H e R and i > 0, 



|<M/«) - 0oo(M)l < I<MaOI + l</>oo(/f)| 



(124) 



and 



< 



:\>K |1 + y]/x]>K V K 



-73 = ^?- ( 125 ) 



On the other hand, notice that (|12ip implies that for r > C 2 + 1, 

|2 + (r-l)(2-0 o (x))| >e 3 r (126) 

for some e 3 depending only on C 2 and e 2 . Then, using the explicit 
expression (TT08j) . eq. (fT26|) and eq. (fT20l) . 



, AW I<1^<£> (T > C2 + 1) , 

£3 t e 3 |^| 



(127) 
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so, for r > C 2 + 1, 



n\>K 



C 2 

\(j> t (n)\ 2 dn < 16-i 

e 3 ^|/i|>-K 



-^=<32- 



C? 1 



(128) 



Equations (|125p and (|128f) together now give 

./|u|>tf ^ V ^3 / 



(r>C 2 + l, A">0). 

(129) 



Step 4: Estimate for |/i| small. Denoting x = fi/r, the differ- 
ence 4>t(fJ>) — 4>oo{v) can be written as 



0t(M) - </>oo(a0 



2i/x(0o(x) - 2) + 2/i 2i + 



<^o(x-)-2 



(l + iM)(2 + (l-T)(^,(x)-2)) 
For the numerator we use (|113|) and (|116l) to get 

2i/x(0 o (x) - 2) + 2/i ( 2i + ~ 2 



(130) 



< (4 + M 2 ) |x| = (4 + M 2 ) | M | 2 ~ (131) 

r 

For the denominator, we use eq. () 1 1 9[> : 

|2 + (l-r)(A)(x)-2)|>i(l + H) (|^| < e a ). 



Then, for |^| < ei, 



|^(m) - 0oo(m)I < C3 



1 H' 



<c 3 -, 



'r(l + |H) 2 - V 
with C3 := 4(4 + M 2 ). Hence for any K := e±r we have 



(132) 



1 



1 



2C 3 2 ^ = 2C 3 2 e 2 -i=:C4^. (133) 



Step 5: Final estimate. From ([129]) and (fl33l) . taking K := eir 
we obtain for r > C 2 + 1: 

/+00 1 / C 2 \ 

- 0oo(^)| 2 < ^ [Ci + 2 + 32-^-J =: Cge"*. (134) 

This shows (|109|) . and finishes the proof. 

□ 
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